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Flight  Control  of  Hypersonic  Scramjet  Vehicles  Using  a 
Differential  Algebraic  Approach 

Tony  A.  Adami*  and  J.  Jim  Zhu  ^ 
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Michael  A.  Bolenderf ,  David  B.  Doman§ ,  and  Michael  Oppenheimer^ 

Air  Force  Research  Laboratory,  Wright  Patterson  AFB,  OH  45433 

[Abstract]  Trajectory  Linearization  Control  (TLC)  is  applied  to  a  longitudinal 
hypersonic  scramjet  vehicle  (HSV)  model.  The  TLC  algorithm  is  based  on  Differential 
Algebraic  Spectral  Theory  (DAST)  which  features  a  time-varying  eigenvalue  concept  and 
avoids  the  use  of  so-called  frozen-time  eigenvalues  that  can  lead  to  unreliable  results  when 
applied  to  time-varying  dynamical  systems.  A  TLC  controller  was  first  designed  for  a 
nonlinear,  affine,  rigid-body  model  using  an  allocation  strategy  based  on  trim-condition 
lookup  tables.  The  tables  were  populated  by  trimming  the  model  at  multiple  operating 
points  while  varying  velocity  and  altitude.  The  trim  data  was  then  fitted  to  a  cubic 
polynomial  function,  and  the  lookup  tables  were  replaced  by  analytical  expressions  for  the 
effector  settings.  The  TLC  design  was  then  verified  on  a  first-principles  based,  longitudinal, 
rigid-body  hypersonic  vehicle  model,  and  initial  simulation  testing  results  are  presented. 


Nomenclature 


X 

=  State  vector 

n 

=  Control  vector 

V 

=  Output  vector 

h 

=  Altitude 

7 

=  Flight  path  angle 

vt 

=  Velocity 

e 

=  Pitch  angle 

Q 

=  Body  pitch  rate 

a 

=  Angle  of  attack 

L 

=  Total  lift 

T 

=  Total  thrust 

My 

=  Pitch  Moment 

-^aero 

=  Aerodynamic  lift 

-^aero 

=  Aerodynamic  drag 

1  prop 

=  Propulsion  force 

-^prop 

=  Lift  due  to  propulsion 

-^prop 

=  Drag  due  to  propulsion 

-^ctrl 

=  Lift  due  to  elevator  deflection 

-^ctrl 

=  Drag  due  to  elevator  deflection 

M0 

=  Mach  number 

<5th 

=  fuel  equivalence  ratio 

Se 

=  elevator  deflection 

PhD  Candidate,  Electrical  Engineering  and  Computer  Science,  Ohio  University,  Student  Member  A1AA 
f  Professor,  Electrical  Engineering  and  Computer  Science,  Ohio  University,  Professional  Member  A1AA 
|  Aerospace  Engineer,  Senior  Member  A1AA 
§  Senior  Aerospace  Engineer,  Associate  Fellow  A1AA. 

Electronics  Engineer,  Member  A1AA 


1 

American  Institute  of  Aeronautics  and  Astronautics 


I.  Introduction 


Hypersonic  Scramjet  Vehicle  (HSV)  technology  has  become  attractive  in  recent  years  as  a  possible  launch  vehicle 
solution.  Current  launch  costs  (approximately  $  1 0,000/lb)  are  prohibitive,  and  Scramjet  vehicles  offer  significant 
relief  by  utilizing  oxygen  in  the  atmosphere  for  combustion,  thus  greatly  reducing  launch  weight.  A  combined-cycle 
approach  is  eventually  envisioned  using  a  vehicle  that  can  alter  its  propulsion  method  based  on  flight  conditions. 
High-speed  travel  and  delivery  of  cargo  (i.e.  global  overnight  delivery)  is  another  potential  application  for  HSV's. 

Research  and  development  in  this  area  is  currently  focusing  on  several  key  challenges  arising  from  the  unique 
integration  of  airframe,  propulsion  system,  and  flight  control  system  of  HSV's.  The  supersonic  velocities  attained  by 
the  air  flowing  through  the  combustor  dictate  a  long,  slender  vehicle  fuselage,  thus  inevitably  leading  to  structural 
flexibility.  The  dependence  of  the  propulsive  efficiency  on  the  geometry  of  the  airframe  for  the  ram  compression 
and  external  combustion  and  expansion  leads  to  a  coupling  between  the  aerodynamics,  structural  dynamics  and  the 
propulsion  system  which  must  be  addressed  by  the  flight  control  system. 

An  analytical  model  based  on  Newtonian  Impact  Theory1'3  and  aero/thermo  analysis  of  the  flow  in  a  Scramjet- 
type  propulsion  system  has  been  among  the  most  commonly  used  for  simulation  and  control  design.  The  dynamics 
are  strongly  influenced  by  both  propulsive  and  aerodynamic  effects,  and  structural  dynamics  of  the  (assumed) 
elastic  vehicle  are  accounted  for  using  a  lump-mass  model.  Using  Lagrangian  methods4,  a  model  has  been  obtained 
which  captures  the  dynamics  due  to  the  rigid-body,  elastic  deformation,  fluid  flow,  rotating  machinery,  wind,  and  a 
spherical  rotating  Earth.  Force,  moment,  and  elastic-deformation  equations  are  derived,  as  well  as  appropriate 
kinematic  equations.  A  model  based  on  Compressible  Flow  Theory  is  under  development  3  that  captures  the  flexible 
structure,  aerodynamic,  and  propulsion  system  coupling.  This  model  will  be  described  below,  and  will  be  used  in 
this  research  for  control  design  algorithm  development  and  testing. 

Study  of  the  flight  dynamics  of  HSV's  has  shown  that  the  vehicles  are  nonlinear,  time -varying,  coupled, 
unstable,  and  non-minimum  phase.  Furthermore,  initial  simulations  indicated  very  slow  response  to  a  reference 
command.  The  structural  dynamics  introduce  further  complications  to  the  control  problem,  including  modal 
uncertainty  due  to  coupling  effects,  potentially  time-varying  modal  frequencies,  and  the  presence  of  structural 
modes  with  frequencies  within  the  control  bandwidth.  Propulsive  effects  can  be  essentially  viewed  as  disturbances 
to  the  flight  control  system.  Note  that  elevator  effects  constitute  an  unstable  zero  dynamics  that  renders  the  overall 
system  non-minimum  phase.*  If  these  effects  can  be  canceled  using  throttle  (possible  if  it  produces  significant  direct 
lift  at  or  forward  of  the  CG  and  responds  quickly)  or  by  using  additional  control  effectors  such  as  a  canard,  then  the 
overall  non-minimum  phase  behavior  will  be  reduced  or  eliminated. 

Flight  controller  design  methods  for  hypersonic  and  other  elastic  vehicles  have  been  an  important  topic,  and 
interest  continues  to  grow.  A  robust  control  system  design  has  been  presented'1'7  that  uses  a  genetic  algorithm  to 
search  a  design  coefficient  space.  Each  search  point  is  evaluated  using  a  Monte  Carlo  routine  that  estimates  stability 
and  performance  robustness.  Using  a  0-1  scale  for  probability  of  the  vehicle  going  unstable,  the  resulting  controller 
reduced  the  likelihood  of  instability  from  0.816  to  0.014  compared  to  open  loop  operation.  Controllers  based  on  the 
mixed-sensitivity  Hoo  technique8  have  been  presented  that  achieve  good  performance  over  a  large  flight  envelope. 
The  authors  address  concerns  of  instability  inherent  in  the  technique  by  ensuring  that  the  frozen  operation  design 
points  are  chosen  carefully.  An  output  feedback  technique9  using  a  novel  error  observer  has  shown  good 
performance  in  the  presence  of  flexible  modes.  An  examination  of  the  flight  dynamics  of  hypersonic  air  breathing 
vehicles10  reveals  the  extensive  coupling  between  flight-path  and  attitude  dynamics  in  the  hypersonic  region. 
Additionally,  this  study  demonstrates  that  propulsive  performance  varies  with  altitude,  Mach  number,  and  angle  of 
attack.  The  presence  of  non-minimum  phase  characteristics  of  the  lifting  body  design  are  addressed11  by  proposing 
an  additional  control  surface.  The  presence  of  a  canard  that  is  ganged  with  the  elevator  can  force  the  right-half  plane 
transmission  zero  further  to  the  right,  allowing  for  increased  freedom  in  control  design  in  the  form  of  increased 
controller  bandwidth.  A  pseudo  dynamic  inversion  (DI)  controller  is  proposed12  that  uses  an  additional  feedback 
loop  to  stabilize  the  zero  dynamics.  This  avoids  the  right-half  plane  pole  that  DI  would  introduce  in  an  attempt  to 
cancel  the  non  minimum  phase  zero.  Another  Dl-based  method  uses  a  spline-interpolated  lookup  table13  during 
simulation  to  calculate  the  force  coefficients,  and  has  shown  promising  results. 

In  the  current  research,  the  authors  study  the  application  of  Trajectory  Linearization  Control14  10  (TLC)  to  the 
HSV  control  problem.  Nonlinear  tracking  and  decoupling  control  using  TLC  is  based  on  Differential  Algebraic 
Spectral  Theory  (DAST),  and  can  be  viewed  as  the  ideal  gain-scheduling  controller  designed  at  every  point  on  the 
flight  trajectory.  It  is  an  analytical  design  that  provides  robustness  by  design  (as  opposed  to  worst-case  optimal).  It 
combines  open-loop  and  closed-loop  control,  and  is  model  based,  thereby  making  maximal  utilization  of  known 
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vehicle  properties.  In  addition,  it  is  adaptive,  using  time -varying  bandwidth  (gain)  to  allow  adaptation  to  modeling 
errors  and  uncertainties,  disturbances,  control  saturations,  and  failures. 

The  inherent  flexibility  of  the  airframe  is  an  especially  vital  consideration,  and  research  in  this  area  is  very 
active  as  well.  An  online  adaptive  filtering  technique17  that  can  simultaneously  identify  multiple  modes  for  flexible 
parameters  has  been  introduced,  although  results  are  limited  to  widely-spaced  frequencies.  A  procedure  based  on 
deeper  understanding  of  acuator  dynamics  is  suggested18,  and  results  in  a  less  conservative  design. 

This  paper  describes  the  HSV  models  under  study  in  Section  II,  and  provides  the  theoretical  background  and 
practical  application  of  TLC,  including  the  controller  design  process  for  the  rigid  body  HSV  longitudinal  dynamics 
model  in  Section  III.  A  control  allocation  scheme  is  discussed  in  Section  IV.  Section  V  describes  the  current 
Matlab/SIMULINK  implementation  of  TLC  on  the  AFRL  model,  and  the  trim-testing  and  step  response  results  are 
presented.  Section  VI  concludes  the  paper  with  a  summary  of  the  main  results. 

II.  Modeling 


A.  First-Principles  Based  Hypersonic  Vehicle  Model 

Researchers  at  Wright  Patterson  Air  Force  Base  (WPAFB)  in  Dayton,  OH  are  developing  a  generic,  hypersonic 
vehicle  longitudinal  dynamics  model  that  takes  into  account  structural,  aerodynamic,  and  propulsion  system 
coupling.5  Compressible  flow  theory  is  used  to  determine  the  pressures  acting  on  the  vehicle,  and  oblique  shock 
theory  is  used  to  determine  the  angle  of  the  bow  shock.  The  forces  acting  on  the  airframe  are  derived  in  the  body 
axis  as  Fx  and  Fz,  where  the  propulsion  and  aerodynamic  forces  and  moment  My  are  heavily  coupled  and  are 
nonlinear  functions  of  the  state  variables  altitude  h,  velocity  Vt,  angle  of  attack  a,  and  control  variables  elevator 
deflection  <5e,  engine  throttle  (equivalence  fuel  ratio)  <5th,  and  the  cowl  door  distance  x&.  The  model  has  been 
implemented  in  MATLAB,  and  will  be  treated  as  a  black  box. 

B.  Nonlinear  Affine  Rigid  Body  Model  for  TLC  Design 

Typical  modern  nonlinear  flight  control  techniques  employ  the  affine  nonlinear  state  equation  model 

x  =  f(x)+g(x)n 

Tj  =  h,(x)  +  d{x)n 

where  x  is  the  state  vector,  fj,  is  the  control  vector,  and  ??  is  the  output  vector.  The  vector  field  /  usually  captures 
inertial  and  structural  couplings,  as  well  as  aerodynamic  and  propulsion  couplings  of  the  state  variables,  while  the 
vector  field  g  represents  effectiveness  of  the  control  effectors  on  the  rate  of  change  of  the  state  variables. 

In  order  to  maximize  the  crosscutting  capability  of  our  design  for  scaling  or  alteration  of  the  airframe,  or 
migration  to  different  airframes,  we  choose  to  seperate  the  inertial/structural  dynamics  from  the  aero/propulsive 
dynamics  as 

x  =  fi(x)  +  f2(x)  +  g(x)/i 

where  fi(x)  captures  the  inertial/structural  dynamics  due  to  vehicle  mass  properties  and  strucural  configuration  that 
can  be  determined  by  analysis  and  testing,  while  ^(x)  represents  aerodynamic  and  propulsive  forces  and  moments 
that  are  usually  determined  by  CFD  analysis,  wind  tunnel  testing,  and  flight  data.  The  control  design  will  first  use 
f 2  ( x )  as  virtual  controls  to  achieve  desired  motion,  then  realize  the  virtual  control  commands  as  a  (static)  control 
allocation  design.  To  this  end,  and  at  the  present  stage,  the  longitudinal  rigid  body  equations  of  motion  in  the 
stability  frame  is  given  by 


h  =  Vt  sin7 


7  = 


in 


g  COS7 


1 

Vt 
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T 

Vt= - g  siiry 

m 

e  =  Q 

Q  =  - -My 
1yy 

where  L  is  the  total  lift  force,  T  is  the  net  thrust  force,  and  My  is  the  total  pitching  moment.  These  variables  will  be 
used  as  virtual  controls  L  =  /ii,  T  —  (i 2,  My  =  and  are  decomposed  into  components 

L  =  -^aero  T  -^prop  T"  -^ctrl 
T  =  —Daero  T*  Tpr0p  Dctr\ 

My  =  My^QYQ  +  -^V/.prop 

The  physical  constant  parameters  are  given  by: 

m  =  300  slug,  g  =  32.17  ft/s2,  i?e  =  2.093  x  107,  Iyy  =  5  x  105slug-ft2. 

With  the  state,  virtual  control  input,  and  output  variables  defined  as  given  in  the  nomenclature,  the  state  and  output 
equations  can  be  written  as: 

h  r  Vt  sin  7  1  0  0  0 

1  II  n  n  r  n 

7  -g  y  COS7  m  vt  u  u  L 

Vt  =  —  gsin7  +  0  ^  0  T 

9  Q  0  0  0  [  My. 

Q  [  0  J  OOf 

=  f{x)+g{x)y, 


L  //prop  “I  f/aero  T"  Tctri 

T  =  Tprop  -  Daero  ~  Dctri  =  b(6,  a,  Vt) 

Vfy  /If, prop  T  J///y.aero 


Note  that  in  this  simplified  case  /i  is  neither  regulated  nor  measured.  However,  the  altitude  will  have  an  effect 
on  the  lift,  drag,  propulsion,  and  control  forces  and  moments  due  to  the  variations  in  density  and  temperature,  and 
will  be  included  in  the  design.  For  autonomous  flight,  the  regulated  variables  can  be  the  altitude  and  down  range 
trajectories.  In  that  case,  flight  path  angle  and  velocity  become  intermediate  variables. 


C.  Nonlinear  Affine  Virtual  Control  and  Control  Effector  Model 

This  study  assumes  an  affine,  nonlinear  actuator  model  for  virutal  control  and  control  effectors  in  the  form  of 

P(X )  =  Ptrim  +  P 


where  ptrim  =  pftrim)  is  the  nominal  value  of  p  at  the  trimmed  (equilibrium)  flight  condition  xtrim,  and  p  is  the 
perturbation  from  the  nominal  which  will  be  approximated  by 


x 

■frim 
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where  x  is  the  perturbation  of  the  flight  condition  from  xtrim-  The  forces  and  pitch  moment  are  given  by 

Taero  =  C La&  H-  -Taero?trim 

T^aero  =  C Da  &  T"  T)aero?trim 

-Tetri  =  CLSe8e  -f  L ctrl, trim 
T^ctrl  =  ^Dse^e  T  -f^ctrl,trim 

Tprop  =  C$th ^th  T  -^prop  trim 

Tprop  =  (-^prop  “b  -^prop,trim)  sinQ;  —  -^propSifr^trim  ~b  Tprop,trimSinQ;trim 
Tprop  =  (-Tprop  T"  -Tprop, trim)  COSO;  —  -^propCOSC^trim  T  T^prop^rimCOS^trim 

My  =  C§e8e  +  Myjr[m 


Thus 


T  ^5£/i^th^i^^trim  T"  Cxa ^  T"  ^Lse^e 


T  —  <5thCOSOJtrim  C]j  CX  CDtJ^e 


where 


My  —  CM6{h8\h  +  Cmq&  +  C'Mse^e 


L  —  -T  -Ttrim  —  -T  -^prop,trim  sinCktrim  T/aer09trim  T/Ctrl,trim 


T  —  T  Ttrim  —  T  -Tpr0p  trim  COSQ^trim  T~  T)aero?trim  T"  T)ctri?trim 


My  —  My  My  trim 


—  Q;  O^trim 
=  ^e,trim 


^th  —  ^th  ^th,trim 

In  this  case,  the  parameters  of  the  model  are  given  by 


CW(F£  ,trim?^trim)  — 


dR 


prop 


cMjyt  ,trim5^trim)  — 


96th 

9M„ 


96, , 


t^trirrb^tri 


tritni^tr 
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@La\Vt,  trim^trim)  — 


dLa, 


da 

dDaero 


da 

dMv 


da 

dL  aer0 


CMa(Vt  .trim-^trim  j  — 

CijeC^i.trimj^trnn)  = 

CD4e  (^.trim^trim)  = 


dSc 

dDaer0 


Cmjc  (  trim;  ^'trirn  )  — 


% 

dM„ 


d8e 


in 


Ltnm  .Lnm 


tO.tnm  .^)nm 


^i.tnm/1  ti 


t/  trim.^tr 


where  the  trim  values  and  derivatives  are  obtained  from  the  AFRL  Model  by  performing  a  two-dimensional 
optimization  at  multiple  operating  points  based  on  velocity  and  altitude. 

ITT.  Trajectory  Linearization  Controller  Design  and  Simulation  Implementation 


A.  Trajectory  Linearization  Control  (TLC) 

TLC  provides  robust  stability  and  performance  without  interpolation  of  controller  gains,  and  eliminates  costly 
controller  redesigns  due  to  minor  airframe  alteration.  Controller  gains  for  the  commanded  torque  are  computed 
symbolically  based  on  the  theoretical  derivations  and  are  applicable  to  any  feasible  trajectories  and  airframes  with 
known  mass  properties.  Therefore  the  design  is  scalable  and  mission  adaptable.  Additional  features  include  the  use 
of  correct  time -varying  eigenvalue  stability  theory,  and  automated  design  using  symbolic  toolbox  script.  There  is  no 
interpolating  for  gain-scheduling  needed,  since  the  method  provides  'ideal'  gain-scheduling  at  every  point  along 
guidance  trajectory. 

Figure  1  shows  a  conceptual  block  diagram  of  the  TLC  method.  A  so-called  pseudo-inverse  is  obtained  for  the 
nonlinear  plant  model  to  calculate  the  open-loop  nominal  control.  The  tracking  error  is  then  driven  toward  zero 
exponentially  using  a  linear  time-varying  (LTV)  tracking  error  regulator. 


Figure  1.  Conceptual  Configuration  Trajectory  Linearization  Control. 
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B.  TLC  Guidance  and  Attitude  Controller  Design 

The  control  design  method  has  been  applied  to  the  affine,  nonlinear,  rigid-body  vehicle  longitudinal  model  with  a 
simplified  linear  control  allocation.  The  design  will  be  carried  out  for  a  slow  translational  (guidance)  loop  and  a  fast 
rotational  (attitude)  loop  based  on  time-scale  separation  (singular  perturbation  principle). 


1.  The  Guidance  (Translational  Motion)  Loop 

The  reference  command  inputs  are  the  commanded  flight  path  angle  7com  and  the  commanded  velocity  V)coin .  The 
(virtual)  control  effector  inputs  to  the  guidance  loop  are  chosen  to  be  the  throttle  <5th  and  the  commanded  angle  of 
attack  acom.  The  outputs  are  the  same  as  the  overall  system.  The  effects  of  the  elevator  deflection  6e  are  treated  as  a 
disturbance  that  is  ignored  in  the  controller  design.  The  guidance  loop  equations  of  motion  are  given  by 


1 - 1 

Tv  ->■ 

i _ i 

= 

~9yt  cost’ 
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i  i 

0  " 
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—  /gi(xgi)  +  ggi(xgi)figi 
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Figure  2.  HSV  TLC  Controller  Configuration 
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i.  Nominal  Guidance  Controller  Design: 


Since  h(x)  is  independent  of  the  input  /zgi,  the  relative  degree  is  greater  than  zero.  Taking  the  time  derivative  of  771 
once  yields 

1  1  1  r 

Vi  =  7  =  ~9  77  COS7  +  —  —  L 
Vt  m  Vt 

>h  =  V  =  :  -5sin7+  —T 
m 

Thus  the  system  has  a  vector  relative  degree  [1  1] .  The  inverse  mapping  from  fj  to  /2gI  is  given  by 

L  =  m  (  Q  7j  1  +  g  COS7 ) 

T  =  m(fj2  +  gsirfy) 

The  first-order  pseudo-differentiators  are  designed  by 

i  =  1,2 

Note  that  rji  =7  ,rj2  =  Vt-  The  state  equation  of  the  pseudo  dynamic  inverse  of  the  plant  is  then  given  by 


7 

_ 

— Wdiff,i72 

+ 

vt. 

Wdiff^I ft. 

Wdiff,l 

0 


—  ^gi(xgi)  +  rgi(xgi)  r ] 


0 

Wdiff,2 


Vi 
V  2 


L 

T 


—m(v  (Wdiff,i7  —  g  COS7  ^ 

+ 

mV  tWdiff.i 

0 

~  Vi~ 

-m[um,2Vt  -  g sitfyj 

0 

fnuj  diff,2 

.  rl2_ 

®gi(®gi)  +  f?gi(a;gi)  77 


Clearly,  the  pseudo-inverse  is  globally  exponentially  stable.  It  is  noted  that  in  unaccelerated  cruise  flight,  the 
nominal  thrust  T  =  0.  This  is  because  the  aerodynamic  drag  force  is  not  included  in  the  plant  model.  Rather,  it  is 
included  in  the  guidance  control  allocation.  Thus,  T  =  0  does  not  mean  the  throttle  <*>th  =  0.  In  fact  =  6th, trim, 
which  generates  the  thrust  that  cancels  out  the  aerodynamic  drag  force. 


ii.  Guidance  Tracking  Error  Stabilization  Controller  Design: 
Define  the  (augmented)  tracking  error  variables  by 


’  7i„t(i)  " 

fo(l(T)-l(T))dT 

72  (t) 

7(f)  -7(f) 

Vint(t) 

fo(vt(T)  -  Q  (T))dT 

.  V(t)  . 

V(t)  —  Q  (f) 

-L(ty 

‘ m-L(ty 

T(t). 
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r*»iW 

to 

1 _ 

V2  (t)-V2(t)_ 

Then  the  tracking  error  dynamics  is  given  by 


"  i'int  " 

- 

72 

~ 

0 

- 9 

[xfecos(7+7(i))] 

+ 

m 

■^_(L  +  I(f))- 11(f) 

K 

vt 

0 

I  ~g  [sin(7+7  (i))  +  sin7(f)]J 

-t2 

m  ^ 

/gi(xgi,  a^gi(t))  +  §gi(^gi,  /igi,  Xg\(t),  Mgi(^)) 


Vi 

’  7  ’ 

Jh. 

A 

35gl(*)) 


Linearization  along  the  nominal  trajectory  yields  the  linearized  tracking  error  dynamics 

Xg\  =  Agi(t)xg\  +  Bg\(t)ug\ 

2/gt  =  CgiW^gi  +  ^gi(^)Mgi 


where 


gl(i)  “  5xgl 


cgi,Mgi 


Bgi(t)  — 


9U 


(9/Xgi 


c«l(f)  -  dig, 


cgi>Mgi 


^gUMgi 


o 


1 


0 


0 


0  5nksin7(i)  0  5^-C0S7(t) 


Q(t) 

0  0  0 

0  —gcos*f(t)  0 

0  0 
IX  0 

m  Q(t) 


Q(t) 


0  10  0 
0  0  0  1 


Let  the  desired  closed-guidance-loop  (CGL)  tracking  error  dynamics  be  given  by 

Acsi(t)  = 


0  10  0 

— ain(t)  -ai12(t)  0  0 

0  0  0  1 

0  0  — ai2i(i)  — ai22  (t) 


where 


otflW  =  u2u(t) 
ai  a{t)  =  2(uuu(t) 


i=  1,2 


u>u(t) 


Uii(t) 

Design  a  proportional-integral  (PI)  state  feedback  control  law 

ttgl  —  li  cgl  ( ^X'gl 

^I,cgl(f) 


—  Kp,csl(t) 

where  the  controller  gain  matrix  is  found  from 


X2 
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X\ 
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-Kcgi  (t)  —  [RgiBgi(t)]  i?gi[Agi(i)  ^-gi(f)] 


mQ  (t) 
0 


0 


— ain(i) 


m 


0 


-ain(t)  +  g  ^  sin'y  (t) 
5  cosy  (f) 


0  9Wwcosl,{t) 

—  <3.121  (t)  — «122  (t) 


where 


Re  i 


1  0 
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0 
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iii.  Guidance  Loop  Control  Allocation 

The  overall  guidance  loop  virtual  control  is  given  by 

Mi 
M  2 


in 

U2 


The  guidance  loop  control  effectors  are  chosen  to  be  the  fuel  equivalence  ratio  (throttle)  <5*  and  the  angle  of  attack 
a  (in  degrees). 


B.  The  Attitude  ( Rotational  Motion)  Loop 

The  reference  command  is  the  commanded  angle  of  attack  acom  generated  by  the  guidance  loop.  However,  the 
attitude  loop  state  variables  are  chosen  to  be  the  Euler  pitch  angle  6  and  the  body  pitch  rate  Q,  which  are  also 
measured  for  state  feedback.  Thus,  the  commanded  angle  of  attack  is  achieved  via  the  commanded  pitch  angle 
$Com  =  CTcom  +  7-  The  virtual  control  input  is  chosen  to  be  the  pitch  moment  My  =  My.  The  regulated  output  for  the 
attitude  loop  is  the  pitch  angle  6.  The  equations  of  motion  for  the  attitude  loop  are  given  by 


'  6  ' 

'o' 

+ 

'  0  " 
1 

_Q_ 

0 

.  Iyy  _ 

[My] 


=  /alOal)  +  3al(a:)/Xal 

M  =  0  =  h3i(x  ai) 


i.  Nominal  Attitude  Control  Design 

The  attitude  EOM  is  a  linear  system  with  relative  degree  2.  Thus  the  inverse  mapping  from  tj3  to  Ji:,  is  given  by 

M3  =  ^yyV  3 

The  second-order  pseudo-differentiators  are  designed  by: 

^3  =  wdiff,3Wi  -  (2Cdiff,3Wdiff,3^3  +  ^diff,3^3) 

Note  that  fj3  =  Q ,  r[3  =  9  .  The  state  equation  of  the  pseudo  dynamic  inverse  of  the  plant  is  then  given  by 
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^3  J 


±yy 


2  Cdiff,3  Wdiff,3$  +  +  [  4j/Wdiff,3  0  3 


=  ^(Sal)  +  ^al(Sal)  V 3 

Clearly,  the  pseudo-inverse  is  globally  exponentially  stable. 


ii.  Attitude  Tracking  Error  Stabilization  Controller  Design: 
Define  the  (augmented)  tracking  error  variables  by 


0int  (t) 
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Then  the  tracking  error  dynamics  is  given  by 
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Linearization  along  the  nominal  trajectory  yields  the  linearized  tracking  error  dynamics 

Xa\  =  ^4_al  ( -(-  Ba]  (t)Ua| 

2/al  =  Cai(t)xai  +  C>al(i)wal 
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Let  the  desired  closed-attitude-loop  tracking  error  dynamics  be  given  by 


A;al4)  — 


0 

0 

-a2i(t) 


1  0 

0  1 

-a22(t)  -a23(f) 


where  for  constant  closed-loop  eigenvalues  pi,p2,  p3, 
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«21  (t)  =  —P1P2P3 

«22  (t)  =  PlP2  +  p2p3  +  P3P1 

«22  (i)  =  ~(pi  +  P2+  P3)- 

The  proportional-integral  (PI)  state  feedback  control  law  is  designed  by: 

Ua\  =  ft  cal  (  f  )  -X'  a  I 

where  the  controller  gain  matrix  is  found  from 

A"cal(i)  =  [AaiBal(f)]_1f?ai[2leai(f)  -  Aa\(t)} 

=  [Iyy][-a2i(t)  — a22(f)  -a23(f)] 

where 

Aal  =  [0  0  1] 


iii.  Attitude  Loop  Control  Allocation 
The  overall  attitude  loop  control  is  given  by 

[  My  ]  =  [  My  ]  +  [  tial  ] 

The  allocation  will  be  designed  together  with  the  guidance  loop  in  the  next  section. 

IV.  Control  Allocation  Design 


Inverse  Jacobian  Approach 

The  control  allocation  is  redesigned  by  inverting  the  equations  for  L,  T,  and  M.  Define  the  Jacobian  matrix  as 
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Assuming  the  Jacobian  is  nonsingular,  we  have 
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6th 
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To  facilitate  the  initial  control  allocation  design,  we  use  the  inverse  Jacobian  of  the  following  simplified  Jacobian 
matrix  that  allows  the  control  allocation  to  be  decoupled  into  the  guidance  tracking  loop  and  attitude  control 
tracking  loop,  and  to  circumvent  the  non-minimum  phase  effects  of  the  lift  due  to  the  elevator  that  would  likely  lead 
to  instability. 
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Control  Allocation  Implementation 


Lookup  Tables 

As  a  first  step  to  accommodate  varying  flight  conditions,  a  lookup  table  is  formed  that  contains  trim  values  for  lift, 
thrust,  moment,  and  effector  configurations  based  on  velocity  and  altitude.  This  approach  gives  good  results  and 
will  be  used  a  baseline  performance.  However,  use  of  tables  requires  significant  onboard  data  storage.  For  a  large 
flight  envelope,  this  may  be  cumbersome. 

Curve  Fitting 

An  alternate  method  of  implementing  the  control  allocation  design  above  is  to  perform  curve  fits  to  obtain  analytic 
expressions  for  the  effector  settings.  Using  Gaussian  elimination,  the  data  obtained  for  the  lookup  tables  are  fitted  to 
the  expressions  of  the  form 

6i  =  C\M0  +  c%M2  +  c%Mq  +  C4/1  +  c$h2  +  c§, 

and  the  tables  are  replaced  with  functions  defined  by  these  expressions.  In  the  next  section,  simulation  results  are 
presented  that  utilize  these  fits  for  computation  of  effector  settings.  It  is  noted  that,  while  curve  fitting  methods 
reduce  the  amount  of  onboard  data  storage,  they  do  increase  the  computation  time. 

V,  Simulation  Implementation  and  Results 


AFRL  Rigid-Body  Longitudinal  Vehicle  Model 

In  the  current  work,  the  affine  model  has  been  replaced  by  the  rigid-body  AFRL  model.  The  TLC  design  used  above 
is  applied  directly  and  preliminary  results  using  the  curve  fitted  allocation  parameters  are  presented.  The  controller 
designed  for  affine  model  is  used  to  perform  trim  and  step  response  tests. 

Trim  testing  result 

In  the  initial  simulation,  the  commanded  trajectory  is  trimmed  flight  at  Mach  8.0  and  85,  000  ft.  Figures  4-11  show 
that  the  vehicle  trims  successfully  within  100s. 

Step  response  testing  result 

As  shown  in  figures  12-23,  a  increase  in  velocity  of  500  ft/s  is  commanded  at  f  =  100s.  It  is  seen  that  the  controller 
designed  for  the  affine  model  performs  well  on  the  nonlinear,  rigid-body  vehicle  model. 
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VI.  Conclusion 

Initial  study  of  the  longitudinal  hypersonic  scramjet  vehicle  model  under  development  at  AFRL  has  revealed  a  rich 
but  complex  control  design  problem.  The  tight  integration  of  propulsion  and  airframe  presents  unique  challenges, 
and  makes  addressing  the  structural  dynamics  of  the  vehicle  a  critical  issue.  In  the  paper,  a  TLC  design  is  presented 
and  verified  with  successful  trim  and  step  response  simulation  on  the  rigid-body,  longitudinal  AFRL  HSV  model. 
Current  work  addresses  the  rigid  body  AFRL  model  directly,  and  a  comparison  is  drawn  between  allocation 
schemes  bases  on  lookup  tables  and  fitted  expressions.  Initial  simulation  results  are  presented  and  show  a  stable 
controller  that  can  fly  trim  flights  as  well  as  perform  basic  maneuvers. 
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Figure  3.  MATLAB/Simulink  Implementation  of  the  Rigid-Body  Vehicle  Model  with  TLC  Controller 
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Trim  Tests 

In  the  following  figures,  the  vehicle  is  trimmed  at  Mach  8.0  and  85,000  ft. 

Pitch  Rate  Trim  Testing 
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Figure  4.  Velocity  Trim  Test 
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Figure  7.  Pitch  Rate  Trim  Test 
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Figure  5.  Altitude  Trim  Test 
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Figure  6.  Pitch  Angle  Trim  Test 


Figure  8.  Flight  Path  Angle  Trim  Test 
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Figure  9.Angle-of-Attack  Trim  Test 
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Elevator  deflection  Trim  Testing 
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Figure  10.  Elevator  Deflection  Trim  Test 

Step  Response  Tests 

In  the  following  figures,  the  vehicle  is  trimmed  at  Mach 
command  is  issued  to  increase  the  velocity  by  500  ft/s. 
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Figure  11.  Fuel  Ratio  Trim  Test 
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Figure  12.  Velocity  Response  to  Velocity  Step 
Command 
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Figure  14.  Pitch  Rate  Response  to  Velocity  Step 
Command 
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Figure  15.  Pitch  Angle  Response  to  Velocity  Step 
Command 


Figure  13.  Altitude  Response  to  Velocity  Step 
Command 


18 

American  Institute  of  Aeronautics  and  Astronautics 


x  10  3  Flight  Path  Angle  Response  to  500  ft/s  Step 


Figure  16.  Flight  Path  Angle  Response  to  Velocity 
Step  Command 
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Figure  17.  Angle-of-Attack  Response  to  Velocity 
Step  Command 
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Figure  18.  Elevator  Deflection  Response  to  Velocity 
Step  Command 
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Figure  19.  Fuel  Ratio  Response  to  Velocity  Step 
Command 
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Figure  20.  ^  Response  to  Velocity  Step  Command 
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Figure  21.  ^  Step  Response 
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Figure  22.  Step  Response 
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Figure  23.  Step  Response 
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